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, Abstract 



The self-gravitational correction to a localized spherically symmetric static energy distribution is obtained 
from an upgraded Newtonian model which is energetically self-consistent, and is also obtained from the 
Birkhoff-theorem extension of the unique "Newtonian" form of the free-space Schwarzschild metric into 
fP , the interior region of its self-gravitationally corrected source. The two approaches yield identical results, 

^ ' which include a strict prohibition on the gravitational redshift factor ever being other than finite, real and 

positive. Consequently, the self-gravitationally corrected energy within a sphere of radius r is bounded 
by r times the "Planck force", namely the fourth power of c divided by G. Even in the absence of 
spherical symmetry, the same bound, to within a factor of two, is shown to apply. That energy bound 
rules out any physical singularity at the Schwarzschild radius, and it also cuts off the mass deviation of 
any interacting quantum virtual particle at the Planck mass. Because quantum uncertainty makes the 
minimum energy of a quantum field infinite, such a field's self-gravitationally corrected energy essentially 
attains the Planck force times that field's boundary radius r. Roughly estimating r as c times the age of 
^ . the universe yields a "dark energy" density of 1.7 joules per cubic kilometer. But if r is put to the Planck 

U ' length appropriate to the birth of the universe, that energy density changes to the enormous Planck unit 

value, which could quite conceivably drive primordial "inflation" . The density of "dark energy" decreases 
C/3 . as the universe expands, but more slowly than the density of ordinary matter decreases. Its evolution 

O ' suggests "dark energy" has inhomogeneities, which may be "dark matter". 

r-| ! Self-gravitational correction of spherically symmetric localized static energy 

Ph. a common tacit idealization made in non-gravitational theoretical physics is to ignore self-gravitational 

corrections as being vastly too small to matter. The ratio of the electrostatic to the gravitational attraction 
, of electron to proton is a great many orders of magnitude, for example. But the uncertainty principle 

of the quantum theory can manifest a disconcerting predilection to throw up infinite energies, and if we 
understandably quail at abandoning so firmly established a principle, it behooves us to at least try to ponder 
' its self-gravitational implications. Even in gravity theory itself , well-known solutions for source-free regions, 

, such as that of Schwarzschild, exhibit troubling behavior near an idealized gravitational source which is 

sufficiently strong and compact. That raises the question of whether we have adequately modeled how 
self-gravitational corrections might modify so strong and compact an idealized gravitational source. 
CNJ ' Here we essay a simplified inroad into self-gravitationally correcting localized energies that are obtained 

or postulated under the tacit assumption that G = 0. To make exact calculation readily feasible, we re- 
strict the localized energies actually treated to be spherically symmetric and static. Such models turn out 
to be readily formally solvable in a Newtonian gravitational framework which has been upgraded through 
^ incorporation of the concept that energy (not mass) is to be self-consistently regarded as gravity's source, 

, notwithstanding the complication that gravity unavoidably alters this selfsame energy. The consequences of 

this relatively straightforward energetically self-consistent Newtonian approach are presented in the remain- 
der of this section. In the next section we show that a somewhat less straightforward metric-based general 
relativistic approach yields identical results, albeit that latter approach manifests especially transparently, 
due to its inherent metric nature, the fact that the Schwarzschild singularity isn't physically realizable. 

Starting from a given spherically symmetric static cumulative energy distribution EG=o{r) that is a feature 
of a gravitationally bereft (i.e., G — 0) "world", our goal is to calculate its self-gravitational modification 
Ecir) that is intended to apply to a more "realistic" world where G > 0, but we don't go so far as to permit 
Ec{r) to become nonstatic. In view of its cumulative nature, we assume that EG=o(r) vanishes at r = 0, 
i.e., 

Ea^r = 0) = 0, (la) 
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and that it has nonnegative derivative, 

E'G=o{r) > 0. (lb) 

We shall also implicitly assume that EQ=Q{r) is a physically realistic localized cumulative energy distribution, 
namely that it satisfies E'Q^Q(r) = for all r which are greater than some positive bounding radius R. We 
also note here that E'Q^Q{r) is, of course, closely related to the nonnegative spherically symmetric energy 
density Tc^oir) via, 

E'G=o{-r) = 47rr2TG=o(r). (Ic) 

Now if it were tlic case that we actually had EQ{r) in hand, we could calculate [EQ{r + dr) — EQ{r)], 
namely the gravitationally-modified energy of the immediately enveloping infinitesimally-thick spherical shell, 
by proceeding to discount from that shell's original energy E'Q^f^{r)dr the negative Newtonian gravitational 
work —[GEG{r)/{c^r)]E'Q^Q{r)dr which is required to assemble it at radius r in the presence of the attractive 
gravitational field produced by Ecir). Assembly of that infinitesimally-thick shell at radius r also entails 
negative Newtonian gravitational work done by the shell's constituents on each other, but that work is of 
second order in dr, and therefore can safely be ignored here. Thus we obtain, 

EG{r + dr) - Ecir) = (1 - [GEG{r)/{c*r)])E'Q^o{r)dr. (2a) 

In addition, the fact that EG=o{r = 0) vanishes eliminates any energy carry-through in the degenerate r = 
special case, and leaves Ecir = 0) with no option but to vanish as well, i.e., 

Ecir = 0) = 0. (2b) 

From Eq. (2a) wc see that the physical core of the calculation is the dimensionlcss Newtonian gravita- 
tional energy-reduction factor (1 — [GEG{r)/{c'^r)]) at radius r, which upon multiplication into the radius-r 
infinitesimally-thick enveloping shell's original energy E'Q^Q{r)dr, yields that shell's gravitationally-modified 
energy [Eair + dr) - Ecir)] = E'Q{r)dr. 

In the limit dr 0, Eq. (2a) becomes the linear inhomogeneous first-order differential equation, 

E'a{r) + {Gld'){E'a^,{r)lr)EG{r) = E'a^,{r), (2c) 

for Ecir), with Eq. (2b) as its single boundary condition. Eq. (2c) is readily reduced to quadrature after 

multiplying it through by the integrating factor e'^^'^ ^ Iro^^'^=°^'^ ^^'^ ^"^"^ ^ which, together with Eq. (2b), 
yields, 

P . N (G/c") f iE'^^o(r')/r')dr' .r , , „, . (G/c*) J"' {E'a^^{r")/r")dr" (2d) 

EG(r)e •'^0 = dr EQ^Q{r )e ■'^o , y^"-^) 

a result which is more simply written as, 

£;G(r)=/;dr'£^=o(r-')e-'^/^'^^-'''^^-^'^"^/''"^*^". (2e) 
In light of Eq. (la) we obtain from Eq. (2c) that, 

hm i?G(r) = /; dr' E'^^„{r') = Ea^r), (2f) 

and in light of both Eqs. (lb) and (la), we obtain from Eq. (2e) that, 

< Ecir) < /; dr' E'^^,{r') = EG=o{r). (2g) 

The results given by Eqs. (2f) and (2g) are entirely expected, but in addition to these there lies concealed 
in the deceptively humdrum form for Eci^r) which is given by Eq. (2c) a startling physical feature. To expose 
that property of E^r) to the light of day, we note that the integrand on the right-hand side of Eq. (2e) 
comes remarkably close to being a perfect differential. Therefore Ecir) can be rewritten, 



EG{r) = (cVG) f^dr'r'd (^e-^°''''^ ^r'^'''^^'^^''"^'''"^''''"^ /dr', (3a) 



which neatly lends itself to integration by parts, with the result, 



Eair) = {c^r/G) 
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(l/r)/;dr'e-(^/'=')X'(^--('-")/'-")*" 



(3b) 



In conjunction with Eq. (lb), Eq. (3b) implies that, 

< EG{r) < (cV/G), (3c) 
which can be combined with Eq. (2g) to yield, 

< Ecir) < min [(cV/G), EG=o{r)] ■ (3d) 

We therefore have the fascinating result that regardless of how large the cumulative energy -Bg=o ('') rnay be 
before it is self-gravitationally corrected, its self- gravitational modification Ecir) cannot attain the "Planck 
force" (d^/G) times the radius r of the sphere which encloses it. With this we are at long last in possession 
of a potent instrument with which to confront the pervasive energy infinities that have for so long been the 
quantum theory's dispiriting affront to physical understanding. Nor does the "Schwarzschild singularity" 
seem so troubling in the soothing light of Eq. (3c), which makes it very clear that arbitrarily strong and 
compact gravitational sources simply do not exist. 

Although the "Planck force times radius" localized-energy upper bound of Eq. (3c) may seem unfamiliar, 
it is straightforward to recast it into an equivalent form to which which we can readily relate, namely, 

1 > (1 - [GEG{r)/{c^r)]) > 0, (3e) 

i.e., the dimensionless Newtonian gravitational energy-reduction factor that appears in Eq. (2a) is always 
real and positive. Also equivalently, the inverse of the Newtonian gravitational energy-reduction factor, 
namely the dimensionless Newtonian gravitational redshift factor, is always finite, real and positive. Those 
statements in fact transcend this particular calculation and are endowed with universal validity, having been 
repeatedly robustly verified for the general relativistic gravitational redshift factor under the most extreme 
conditions of gravitational collapse [1]. It is further worthwhile to note here that in the next-to-last two 
sections of this article it will be shown, by using an interesting nonperturbational iteration method, that the 
"Planck force times radius" upper bound on the self-gravitationally corrected energy within any spherical 
region still holds to within a factor of two even in the absence of spherical or any other symmetry for the static 
energy distribution. Therefore in the next few sections wc will apply Eqs. (3c) and (3b) in a broad-brush 
way to some of the self-gravitationally related issues that naturally arise in gravity theory itself, quantum 
field theory and even cosmology. 

Before doing so, however, wc close this section with the technical development of the asymptotic expansion 
oi Ecir) when EG=o{r) greatly exceeds the "Planck force times radius" upper bound {c^r/G) for Eair); that 
asymptotic expansion will serve as a useful reminder of the basic logic which needs to be applied in dealing 
with the energy infinities that can result from quantum theory considerations. 

To develop this asymptotic expansion of Eair) from Eq. (3b), it is clearly useful to define, 

«(r) {G/c'')E'a^oir)/r = 47r(G/c4)rTG=o(r), (4a) 
where the second equality follows from Eq. (Ic). Using K{r), Eq. (3b) reads, 

1 - (1/r) /; dr' e- X' . (4b) 



Eair) = (c^r/G) 



Since K{r") > 0, the exponential integrand of the integral over the variable r' obviously has its maximum 
value of unity at r' = r. Our basic strategy will therefore be to expand the argument of that exponential 
(which itself is an integral over the variable r") around that point r' = r. To facilitate this, we change 
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variable from r' to the dimensionless variable u = {r — r')K{r), after which we can expand the argument of 
that exponential in powers of u. The change of integration variable from r' to u causes Eq. (4b) to read, 



Ecir) = (cV/G) 



K{T")dr" 



(4c) 



Expanding the argument of the exponential in powers of u formaUy produces, 

jr-(u/.(r)) ^ _^ ^ ^(-l)^(u/«(r))^At(^-i)(r)/j! . (4d) 

Putting this formal expansion result into the exponential produces, 

and we now proceed to expand the second exponential on the right-hand side of Eq. (4e), 

eE~.(-i)^W«M)^«"-'M/.-i ^ 1 ^ g (Er=2(-l)^(V«(r))^V^-i)(r)/j!)" /n! . (4f) 

Tl=l 

We now intend to substitute the right-hand side of Eq. (4f) into Eq. (4e), and then to substitute the 
right-hand side of Eq. (4e) into Eq. (4c) . At this point we recall that we are interested in the asymptotic form 
of Eq. (4c) when EG=o{r) > (c^r/G). That state of affairs can be expected to be coincident with the upper 
integration limit rK{r) of the dimensionless u-integration in Eq. (4c) being very much greater than unity. 
Therefore, together with the substitutions mentioned in the first sentence of this paragraph, we shall as well 
change the upper limit of the dimensionless u- integration in Eq. (4c) to oo. That change, together with the 
expansions given by Eqs. (4f) and (4e), enables us to analytically evaluate Eq. (4c) on a term-by-term basis 
by simply applying the elementary formula, 

du e~" vJ' = k\ . (4g) 
Therefore, the first few terms of the desired asymptotic expansion of EQ{r) when rK{r) ^ 1 are given by, 

Ecir) ~ (cV/G) [1 - (l/(m(r))) (l + K'{r)/{K{r)f - K"{r)l{K{r)f + ^{k' {r)f / {K{r)Y + •••)]• (4h) 
In the case that the uncorrected initial energy density Tg=o(^) happens to be constant in r, we have that, 

K'{r)/iK{r)f = l/ir^ir)) = l/i47r(G/c^yTG=o), (4i) 

and also that K"{r) = 0. Therefore in the case of constant uncorrected initial energy density Tg=o we have 
that when r«;(r) » 1, 

Ecir) - (cVG) [1 - (l/(r«(r))) - (l/(r«(r)))2 - 3(l/(r«(r)))3 + •.•], (4j) 

where, 

r«(r) = 47r(G/c4)r2TG=o = 3(G/(cV))£;G=o(r). (4k) 



Self-gravitational correction from a general relativistic metric 

We have strongly emphasized that the physical core of the energetically self-consistent Newtonian self- 
gravitational correction of spherically-symmetric static cumulative energy distributions lies with the dimen- 
sionless Newtonian gravitational energy-reduction factor (1 — [GEG{r)/{c'^r)]) of Eqs. (2a) and (3e). Now a 
metric, which is in effect a system's general relativistic gravitational potential, also provides a dimensionless 
gravitational energy-reduction factor. For a spherically-symmetric static system that metric gravitational 
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energy-reduction factor is (.900 (''))^, which is the inverse of the metric's gravitational time- dilation or red- 
shift factor (500 (f"))"' [2]. Such a spherically-symmetric static system is of course itself described by a 
self-gravitationally corrected cumulative energy distribution, 

EG{r)=AT:j;{r'fTG{r')dr', (5a) 

where Toir) is its corresponding sclf-gravitationally corrected spherically-symmetric static energy density. 
Now Tcir) would be expected to be the energy density of a localized system, i.e., Tair) = for all r > R, 
where R is that localized system's bounding radius. Thus if we denote Ec{R) as simply Eq, we see that 
Eg{i^) = Eg for all r > R, and, in the free-space region, namely r > R, this system's m,etric will be a free- 
space Schwarzschild metric that has the self-gravitationally corrected energy constant Eq. One well-known 
representation of such a free-space Schwarzschild metric with self-gravitationally corrected energy constant 
Eg is its "isotropic" form [3], 

W = ( i;|!g£acV)I )'(^'^^)' - {^HlOEaKc^p)]r{{dpf + {pdef + (psin^^)^). (5b) 

Unfortunately, the gravitational energy-reduction factor (500 (p))^ of this Eq. (5b) "isotropic" form of the 
Schwarzschild metric with energy constant Eg obviously fails to match the corresponding Newtonian grav- 
itational energy-reduction factor (1 — [GEg/{c^p)\). It tmns out, however, that the simple change of radial 

coordinate variable r =^ p -|- 5 [GEg /c^] in the Eq. (5b) "isotropic" form of the free-space Schwarzschild 
metric produces the unique desired form of the free-space Schwarzschild metric, namely, 

[cArf = {i-[GEa/{c^r)]Y{cdtf - {i-^[GEa/(c'r)])-\drf - {i-^[GEa / {c\)])-\{rd9f + {rsm9d<Pf), (5c) 

that actually satisfies {goo{r))^ = (1 ~ ['^^o/ic'^r)]). It is therefore appropriate to refer to the Eq. (5c) form 
of the free-space Schwarzschild metric as its "Newtonian" form. 

Next we wish to extend this "Newtonian" form of the free-space Schwarzschild metric into the interior 
region r < R of its source, where it is no longer necessarily the case that EG{r) is equal to the constant 
Eg- The general relativistic instrument of choice for accomplishing this extension is the combination of the 
Birkhoff theorem and its corollary for spherically-symmetric static systems [4]. 

If we focus on a particular radius value ro, the Birkhoff theorem's corollary tells us that the part of 
the spherically-symmetric static cumulative energy distribution EG{r) for which r > ro makes no contri- 
bution to the local ro value of the metric, while the Birkhoff theorem itself tells us that the part of the 
spherically-symmetric static self-gravitationally corrected cumulative energy distribution EG{r) for which 
r < ro produces at ro a local Schwarzschild metric whose energy constant Eg has the value EG{ro)- 

The upshot of the Birkhoff theorem and its corollary is thus that the "Newtonian" free-space Schwarzschild 
metric form of Eq. (5c) can be extended to the interior region of its self-gravitationally corrected spherically- 
symmetric static source, which is described by the cumulative energy distribution EG{r), by the very simple 
expedient of replacing the energy constant Eg which appears in Eq. (5c) by that cumulative cmergy distribu- 
tion Ecir). This yields the complete "Newtonian" metric for the self-gravitationally corrected spherically- 
symmetric static cumulative energy distribution EG{r), 

{cdrf = {i-[GEG{r)/ic^r)])^{cdtf - {i-i[GEG(r)/ic^r)])-'{drf - {i-L[GEG(r)/{c'r)])-^ {{rdOf + {r sinOdcpf) , 

(5d) 

which exactly corresponds to the dimensionless Newtonian gravitational energy-reduction factor in Eq. (2a) 

{goo{r))i = {l-[GEG{r)/{c^r)]). (5e) 

Therefore if we multiply the gravitationally-modified energy E'Q{r)dr of an infinitesimally-thick spherical shell 

at radius r by the inverse 1/(1 — \GEG{r) / {c'^r)]) of the corresponding radius-r dimensionless gravitational 
energy-reduction factor given by Eq. (5e), we undo the gravitational modification of the energy of that 
infinitesimally-thick shell and arrive at what its energy would be if gravity were "switched off", i.e., at what 
its energy would be if G were put to zero, 

E'G^,{r)dr = E'G{r)dr/{l - [GEG{r)/{c^r)]). (5f) 
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Eq. (5f) implies the linear inhomogeneous first-order diff'erential equation, 

E'air) + {G/c^){E'^^,{r)/r)Ea{r) = E'^^^r), (5g) 

which is exactly the same as the key Eq. (2c) result of the energetically self-consistent Newtonian approach 
set out in the previous section. This differential equation is, of course, completely solved and exhaustively 
discussed in that section, and is shown there to imply the crucial upper hound on the self-gravitationally 
corrected cumulative energy Ec{r) that is given by Eq. (3c), namely that, 

< Ecir) < {c^r/G). (5h) 

The Eq. (5h) upper boimd on the sclf-gravitationally corrected cumulative energy Eoir) is easily verified to 
imply that, despite one's cursory initial impression, the complete "Newtonian" metric given by Eq. (5d) has 
no singularities whatsoever. 

Of course the complete "Newtonian" metric of Eq. (5d) is entirely consistent with the Eq. (5c) "New- 
tonian" form of the free-space Schwarzschild metric that only applies in the r > R free-space region where 
the cumulative energy distribution Ecir) is equal to the energy constant Eg- Therefore once again despite 
one's cursory initial impression, the Eq. (5c) "Newtonian" form of the free-space Schwarzschild metric has 
no singularities whatsoever in its free-space region of validity. Therefore the ostensibly "threatening" sin- 
gularity that occurs in the Eq. (5c) free-space Schwarzschild metric when r equals the Schwarzschild radius 

rs '= GEg/c^ [5] is never physically realized, and this is also the case for the ostensibly pernicious singularity 
that occurs in the Eq. (5c) free-space Schwarzschild metric when r = ^rs. 



Interacting quantum pcirticles that have forbidden energy or mass 

Because of its wave character, a quantum particle which interacts with a potential can penetrate a short 
distance into a region where the potential's value exceeds the particle's energy —such penetration is forbid- 
den to classical particles. In such a region the quantum particle's kinetic energy and momentum squared 
eff'ectively assume negative values, and its penetration length A into that energetically forbidden region can 
be roughly described as, 

A«?i/(-Pe'ff)'- (6a) 

Since for a nonrelativistic interacting particle, 

E = pV(2m) + V, (6b) 

we can rewrite Eq. (6a) as, 

\^h/{2rn{V ~ E))h, (6c) 

and the energetically forbidden region corresponds toV > E. \iV should vary significantly from the edge of 
the forbidden region to the depth A, Eq. (6c) will need to be regarded as an approximate implicit relationship 
which actually needs to be solved for A. That implicit relationship is set up by reexpressing V as a function 
of the distance from the edge of the forbidden region, and the resulting new independent variable is then 
identified as A. 

A highly energetic interacting quantum particle can similarly deviate from the natural rest mass which 
it has when it is free, and thus enter a region of forbidden rest mass. Just as there is an effective length 
limit for quantum particle penetration into a region of forbidden energy, so there is an effective proper time 
limit, i.e., lifetime, for quantum particle penetration into a region of forbidden rest mass, namely, 

T « n/(A„c2), (7a) 

where is the rest-mass deviation experienced by an energetic interacting quantum particle and r is the 
lifetime of that deviant-mass state. Given its limited lifetime, such a deviant-mass virtual particle is as well 
limited in space to a spherical region of radius R k, ct, 

R^CTK h/{Amc). (7b) 
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Given that radius R of Eq. (7b), Eq. (3c) then yields a self-gravitational approximate upper bound on the 
total energy E of this deviant-mass particle, 

E < {c'^R/G) « c^n/{GAm). (7c) 

Since A^c^ cannot exceed the deviant-mass particle's total energy E, it follows from Eq. (7c) that, 

AmC^ < c^h/{GAm). (7d) 

Therefore, 

(A„)2 < hc/G, (7e) 

which implies that, 

Am<{ric/G)i, (7f) 

namely that is approximately bounded by the Planck mass {hc/G)i . 

This approximate Planck-mass upper bound on the mass deviation of any interacting quantum virtual 
particle universally cuts off the ultraviolet divergences which bedevil quantum particle scattering amplitude 
calculations. 

Self-gravitational correction of the infinite energies of quantum fields 

A noninteracting field (e.g., the source- free electromagnetic field) always decomposes into an infinite number 
of independent simple harmonic oscillators whose frequency spectrum has no upper bound. Upon quan- 
tization, each such oscillator has a minimum positive energy (i.e., quantum ground state energy) which is 
equal to its frequency times h/2: that minimum energy is quantum theoretically inviolable, being completely 
mandated by the quantum uncertainty principle. The fact that there are an infinite number of such simple 
harmonic oscillators with an unbounded frequency spectrum implies that the corresponding quantum field 
always has infinite energy. If the field is set up in a bounded region, then its quantum counterpart has not 
only infinite energy, but necessarily infinite average energy density as well. 

Therefore for quantum fields the uncertainty principle baldly confronts us with an unphysical nightmare, 
yet without this selfsame uncertainty principle likewise mandating a definite minimum energy, Rutherford's 
nuclear atom can't be sustained. 

The last apparent hope for the beleaguered theorist in this harrowing circumstance resides with the self- 
gravitational modification of an initially infinite energy density, e.g., Eqs. (4]) and (4k) in the limit that 
the initial energy density Tc=o is infinite. A saving grace of this vertiginous maneuver is that its result is 
particularly unambiguous and simple, namely, 

Ecir) = {cSjG), (8a) 

where r is the bounding radius of the quantum field. Of course it makes no difference whatsoever how 
many or what types of quantum fields reside within that bounding radius: any infinite initial energy density 
produces the result of Eq. (8a). 

The total self-gravitationally modified energy of quantum fields, which is what Eq. (8a) is supposed to 
describe, isn't directly measurable. However, by making use of the fields' bounding radius r, we can obtain 
from Eq. (8a) their self-gravitationally modified averaged energy density p, 

p= (3cV(47rGr2)). (8b) 

It now remains to puzzle out what conceivable "real world" physics could relate to the mathematical ab- 
straction of the quantum fields' bounding radius r. It is, of course, apparent that no material substance 
can serve to bound the arbitrarily high frequencies which such fields are able to muster. The universe' cos- 
mological redshift, however, in principle ought to defang any frequency, and indeed appears to serve as the 
"containment" for all that we can hope to survey. Therefore a not altogether implausible crude estimate of 
the quantum fields' "bounding radius" r ought to be given by the age of the universe [6] times the speed of 
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light, which comes to about 1.3 x 10^*' meters. Putting that vahu; of r into Eq. (8b) yields about 1.7 x 10~^ 
joules per cubic meter (i.e., 1.7 x 10~^ ergs per cubic centimeter or 1.7 joules per cubic kilometer) as a crude 
estimate of the universe' average "dark energy" density. This is in fact of the same order of magnitude as 
what is yielded by observations [7, 8]. 

In addition to the ability of Eq. (8b) to yield a passable crude estimate of the current universe' average 
"dark energy" density, its systematics also seem fascinating. If we project it all the way back to the universe' 
birth, when r was presumbably of the order of magnitude of the Planck length (Gh/c^)^ ^ then p approaches 
of order unity in Planck units of energy density, which is roughly 120 orders of magnitude greater than its 
value for the current universe. 

Tht'orists who did not attempt to actually model the physics which produces self- gravitational energy 
correction have favored this particular enormous value of "dark energy" density because of their adoption 
of a physically-blinkered "universal fix" for infinite results, namely the replacement of any such infinity by 
one Planck unit of the appropriate dimensions [8]. Neither physical modeling of self-gravitational energy 
correction nor observations have much overlap with such undiscriminating replacement of the quantum 
energy-density infinity by its Planck-unit value, but it is still fascinating to consider that enormous Planck 
unit of ''dark energy" density as bciing relevant to the early universe, as that would apparently provide an 
automatic mechanism for the heretofore puzzling "inflation" of that early universe. 

Finally, Eq. (8b) suggests that the average "dark energy" density ought to decrease toward zero as the 
universe; continues its expansion. This brings to mind the not infrequently expressed theorist preference for 
exactly vanishing "dark energy" density over its observed value [8] , which while immensely smaller than the 
Planck unit of energy density, nonetheless absolutely fails to vanish. In fact, completely to the contrary, it 
dominates the net average energy density of our universe [7] . Wc sec that Eq. (8b) apparently caters for all 
tastes in average "dark energy" density, whether those tastes gravitate toward the enormous Planck unit of 
energy density, zero energy density, or anything in between, including a passable rendition of the observed 
average "dark energy" density which actually obtains at the current stage of evolution of our universe. It is 
to be cautioned, however, that while Eq. (8b) indeed has average "dark energy" density decreasing toward 
zero as the universe continues to expand, the average density of normal matter would be expected to decrease 
at a faster rate, so that "dark energy" relative dominance would continue to grow. 

Must self-gravitation be quantized to correct quantum energy infinities? 

Gravity, like electromagnetism, is a gauge theory, and the issues surrounding its quantization formally parallel 
those issues in electromagnetism. In both cases there are dynamical, nondynamical and redundant fields 
present, and the dynamical fields in both cases are two in number and describe transverse radiation. Only 
these two dynamical radiation fields are subject to quantization. 

What remains after the two transverse dynamical radiation fields are accounted for splits evenly into 
nondynamical and redundant fields, neither of which, of course, are subject to quantization. The four- 
potential of electromagnetism yields one redundant field and one nondynamical field of Coulombic character. 
The symmetric metric tensor of gravity yeilds four redundant fields and four nondynamical ones. One of the 
nondynamical fields very roughly corresponds to Newtonian gravity with roughly an energy-density source, 
while the other three merely round out a relativistic four- vector representation, and therefore have roughly 
a momentum-flux source. 

If wc look back at the previous parts of this article, it is clear that the self-gravitational corrections which 
are of overarching importance can all be profitably pondered in a quasi-static or outright static framework. 
The basic ingredients for self-gravitational corrections tend to be Newtonian, albeit an energetically self- 
consistent form of gravitational Newtonianism. 

Gravitational radiation doesn't physically enter into self-gravitational correction, so gravity quantization 
cannot be an issue in such correction, any more than electromagnetic quantization can be an issue in elec- 
trostatics. Gravitostatics is merely more subtle than electrostatics because of its energetic self-consistency. 

Self-gravitational correction of nonsymmetric static energy densities 

In the first two sections of this article we showed that any spherically-symmetric static cumulative energy 
distribution Ec^oir) that satisfies Ea=o{r = 0) = and d{Ec=o{r)) /dr > has a corresponding self- 
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gravitationally corrected spherically-symmetric static cumulative energy distribution Eair) which satisfies 
the inequalities < Ea{r) < EG=o{r) and which also, irrespective of how large EG=o{r) may be, is bounded 
above by the product of its radius r and the "Planck force" (c'^/G) as given by Eq. (3c), i.e., Eair) < r{c^/G). 

We now extend the self-gravitational energy-correction process set out in those first two sections to any 
specified static energy density Ta^oijc) which is nonnegative, smooth and globally integrable, i.e., 

TG=o(r) > 0, (9a) 

Vr (Tg=o{y)) is continuous, (9b) 

and, 

/TG=o(r')dV<oo. (9c) 

Here we specifically refrain from making the assumption of the first two sections of this article that Tg=o{y) 

possesses spherical symmetry, nor do we assume that it possesses aiiy other particular symmetry. 

Now if it were the case that we actually had in hand the self- gravitational correction Tg{y) of the specified 
static energy density TG=o(r), we could calculate the negative Newtonian gravitational work done to bring 
the infinitesimal original static energy TG=Q{r)d^r from infinity to its position at r while subject to the static 
gravitational field that is provided by T(3(r), which yields the result —{G/c^) J (i^r'T(3(r')|r— r'|^^T(3^o(r)'i^r. 
However, because the static gravitational interaction inherently occurs between pairs of infinitesimal energies, 
we must take care to avoid double- counting, so we assign only half of this negative gravitational work 
correction to the infinitesimal static energy located at r, and thus obtain, 

TG{r)d^r = [1 - i(G/c4) /d3r'TG(r')|r - r'|-^] TG=o(r)d3r. (10a) 

Prom Eq. (10a) we see that the dimensionless local gravitational energy-reduction factor J^g{^) that satisfies, 

TG{r)d^r = TG{r)TG^oir)d^r, (10b) 

is given by, 

Tcir) = 1 - i(G/c4) / |r- r'|-iTG(r')(iV. (10c) 

This static local gravitational energy- reduction factor J^g (r) is obviously the inverse of the corresponding 
gravitational time-dilation factor, and thus is equal to (500 (r))^ [2]. If we insert the instruction implicit 
in Eq. (10b) into the right-hand side of Eq. (10c), we obtain the following inhomogeneous linear integral 
equation for J^g{^), 

J-G(r) = 1 - i(G/c*) / |r - r'|-'7G=o(r')J-G(r')rf'r'. (lOd) 
In light of Eq. (9c), we can deduce from Eq. (lOd) that, 

lim J-G(r) = 1. (11a) 

|r|— >(x> ^ ' 

Furthermore, since the integral transform kernel —1/ (47r|r — r'|) is the Green's function of the Laplacian 
operator V^, we in addition deduce from Eq. (lOd) that, 

V2j-G(r) = (27rG/c4)TG=o(r)7-G(r), 

which is readily reexpressed as a zero-energy stationary-state nonrelativistic Schrddinger equation [9] for the 
dimensionless wave function Tg{^), namely, 

(-ftV2/(2m) + V^(r)) J-G(r) =0, (lib) 

whose repulsive potential V{r) is defined as, 

V{r) [nn^G/{mc^)] TG=o(r). (He) 
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Because of Eq. (9c) it is clear from Eq. (11c) that, 



lim V{r) = 0, n Id) 

|r|— >oo 

which, in turn, implies that the largc-|r| limit of .7x;(r) that is given by Eq. (11a) is consistent with the 
Eq. (lib) zero-energy Schrodinger equation. 

Having established the connection of the Eq. (lOd) gravitational integral equation to the Eq. (lib) zero- 
energy stationary-state Schrodinger equation, we now furthermore note that for stationary states of positive 
energy E > this Schrodinger equation becomes, 

(-nV2/(2m) + V{r)) (rlVfi) = E{r\ijE). (12a) 

Prom Eq. (lid) we see that as |r| — >■ oo, Eq. (12a) becomes simply, 

{-n'Vl/{2m)) (r\^E) = E{r\^E), (12b) 

whose solutions include all the dimcnsionlcss plane waves e^P ''/'* for which p satisfies |p|^ = 2mE, as well, of 
course, as the linear combinations of these which comprise the the full set of angularly-modulated outgoing 
and ingoing free spherical waves which have this same scalar wave number k = {2mE)'^ /h [9]. Thus we see 
that Eq. (12a) has a massive inherent solution degeneracy. A useful resolution of this solution degeneracy 
can be achieved by reexpressing Eq. (12a) in an inhomogeneous linear form that can only be satisfied by the 
particular permissible |r| — > oo asymptotic behavior which properly accords with the design of a specified 
experiment. 

That idea underlies the Lippmann- Schwinger inhomogeneous modification of Eq. (12a), which forces 
its wave function to behave as a specifically chosen single permitted plane wave e'^ ''/^ plus only outgo- 
ing angularly-modulated spherical waves in the asymptotic region |r| — > oo where Eq. (12a) is adequately 
described by Eq. (12b). If we denote as (rji/'p ) the solution of Eq. (12a) which satisfies this particular per- 
mitted |r| 00 asymptotic behavior, then the inhomogeneous Lippmann-Schwinger equation that uniquely 
describes (rlV'p) is [10], 

(r|V'+) = e'P-'-/'^ - (r|(-n^V2/(2m) - |p|V(2to) - ie)-^V\il^+). (13a) 

From a static gravitational standpoint the relevant feature of the Eq. (13a) inhomogeneous Lippmann- 
Schwinger equation and its wave function {y\iP^) is that, 

(r|V'^)|p=o = -^G(r), (13b) 

as is seen from comparison of the p = case of Eq. (13a) with Eq. (lOd) — to make this comparison one 
must use Eq. (11c) to express V{y) as the appropriate constants times Tg=q{y), and one must also use the 
fact that the integral transform kernel —1/ (47r|r — r'|) is the coordinate-representation inverse (i.e., Green's 
function) of the Hilbert-space "Laplacian" operator = — |p|^/?i^, namely that, 

(r|(V2)-i|r') = -l/(47r|r-r'|). 



Note that the negative imaginary infinitesimal —ie which appears in Eq. (13a) is unnecessary when p = 0, 
which represents a purely static state of affairs that has no distinguishable outgoing versus ingoing spher- 
ical waves. Indeed the massive solution degeneracy of the stationary-state Schrodinger equation given by 
Eq. (12a) collapses when E = Q. 

Given the Eq. (13b) close relationship of the local gravitational energy-reduction factor !Fg{y) to the 
Lippmann-Schwinger wave function (r|^+), it would seem logical to apply well-known general solution 
methods for Lippmann-Schwinger equations to our gravitational Eq. (lOd). Unfortunately, however, the 
only widely- applied fully general solution method for Lippmann-Schwinger equations is perturbational in 
character, and therefore is inherently subject to failure. 
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The struggle to transcend the perturbational Born trap 

If we bring the second term on the right-hand side of the Eq. (13a) Lippmann-Schwinger equation to its 
left-hand side, we obtain, 

(r|^+) + {r\{K -Ep- ie)-^V\i^+) = e*P■"■/^ (14a) 

where K =' (— fi^ V2/(2m)) is the kinetic energy operator and Ep (|pp/(2m)) is tlic kinetic energy 

c-number scalar that corresponds to the c-number momentum vector p. Taking (r|p) =' e^P ''/'', the formal 
solution of Eq. (14a) is, 

(r|V+) = (r|[l + {K-Ep- ie)-'V]-'\p). (14b) 

The only way forward at this point would seem to be expansion of the inverse of the operator in square brack- 
ets in the well-known Born geometric perturbation series, which involves successive powers with alternating 
signs of the particular operator, 

X = {K-Ep-ie)-^V, 

acting on the momentum eigenstate |p) [11]. If the operator X dominates the identity on the momentum 

eigenstate |p), it is not unlikely that the Born geometric series diverges. One might suppose that in such 
instances one could simply recast the Born geometric expansion to be in powers of the inverse of X, since 
one is formally free to choose either the expansion [1-|-X]~^ = 1 — X + X"^ or the expansion [1 -|- X] ~^ = 

X-i[l-FX-i]-i = X-i-(X-i)2-K(X-i)3 . Most unfortunately, however, since X''^ = V''^ {k - E^, 

the operator X^^ vanishes altogether when applied to the momentum eigenstate |p). This luiariticipatcd 
abrupt setback is a stark warning that hidden snares beset Born-style geometric perturbation expansions for 
the Lippmann-Schwinger equation. 

For the p = case of the Lippmann-Schwinger equation that applies to static gravitation, we are, of 
course, particularly interested in arbitrarily large energy densities Tg=o{y), and therefore in arbitrarily strong 
operators V and X. Thus it is clear that we need to entirely abjure Born-style geometric perturbation 
expansion, but how could that conceivably be accomplished in practice! The only possibility is through 
exploration of nonstandard manipulations of the Lippmann-Schwinger equation. 

Returning to Eq. (14a) we now deliberately shun the elegant and natural factorization of the operator [1 + 
X\ on its left-hand side, which can only lead us down the primrose path to the Born geometric perturbation 
series, and instead opt to forcibly factor that side into two mere functions of the vector coordinaie r, the first 
of which is still {r\tp^) because we aim for a result which is at least akin to Eq. (14b), but the second of which 
is repulsively inelegant, being merely the product of ((r|V'p ))~"^ with the original left-hand side of Eq. (14a). 
There is in fact "method" in that gross ugliness, however, because we can now actually arithmetically divide 
the plane wave e'^''^/^ on the right-hand side of Eq. (14a) by that gauche second factor without resorting 
to any kind of perturbation expansion. A very heavy price has been paid in the coin of gross inelegance, 
but the goal of no perturbation expansion whatsoever has been achieved. To be sure, the almost childish 
manipulations just described haven't extracted any final result from Eq. (14a), what they have produced is 
only a basis for refinement through iteration. It is readily seen, however, that the iteration process is devoid 
of perturbational characteristics; it instead resembles a continued fraction. 

The iteration formula we have just extracted from Eq. (14a) is explicitly, 

(r|4"+i)+) = eWf«/[i + ((r|4")+))-i(r|(^ - Ep - ie)-^V\^^^^+)], (14c) 

for n = 0,1,2,..., where, of course, (r|V'p°^^) = e^P'""/^. Therefore {{r\ipp^~^))~^ is obviously well-defined, 
and the form of Eq. (14c) makes it apparent that for n = 1,2,..., ((r|-0p"''^))^^ is well-defined as well. 
That the iteration formula of Eq. (14c) does not have perturbational characteristics, but rather those of a 
continued fraction is also manifest. 

Finally, on inserting p = into Eq. (14c) we obtain the iteration formula for the gravitational energy- 
reduction factor J^G(r), which is, 

4"+^Hr) = 1/[1 + i(G/c4)(4"\r))-i / |r - r'|-^TG=o(r')4"^r')d3r'], (15a) 
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for n = 0, 1, 2, . . ., where, of course, ^'^"^(r) = 1. 

Since rG=o(r) > from Eq. (9a), TG=o(r) is smooth from Eq. (9b), and / TG=o(r')c^^r' < oo from 
Eq. (9c), it is clear from Eq. (15a) that, 

if 1 > J-^"' (r) > 0, then 1 > (r) > 0. (15b) 

Therefore we can conclude that, 

1 > J^cir) > 0. (15c) 
From Eqs. (9a), (10b) and (15c) we can deduce that, 

< TG(r) < TG=o(r). (15d) 

Prom Eq. (10c) and the fact that the gravitational energy-reduction factor J^g{^) satisfies J^g{^) > we can 
deduce that, 

2(cVG)> J \r-r'\-^TG{r')d^r' = J |r" | "^Tg (r + r")rfV', 
where r" =^ (r' — r). Furthermore, we have that. 



/ \r"\-'TGir + v")d?v" > [ |r"|-'rG(r + r")d^r" > [ TG(r + r")d 

J •'|r"|<i? J\t"\<R 

Therefore from the two foregoing lines of displayed integral inequalities we can conclude that, 

(2i?)(cVG)>/ rG(r + r")dV', (15e) 

J\r"\<R 

namely that the self-gravitationally corrected static energy contained in any sphere cannot exceed the di- 
ameter of that sphere times the "Planck force" {d^/G). Note that this result is completely independent of 
any assumption concerning symmetry properties of the energy distribution. Therefore it is likely to be overly 
conservative in practice. A more practical energy upper-bound estimate can be obtained by simply averaging 
the maximum possible and minimum possible values of |r"| that occur when integrating over the sphere of 
radius R, which yields R/2, and therefore the "rough" bound. 



R{c^/G)> [ TG{v + v")d\", 

J\r"\<R 



which is in line with the result of Eq. (3c), for which spherical symmetry was assumed. 

No doubt the most interesting aspect of what has just been presented here is the unfastening of the 
shackles of the perturbational Born-expansion paradigm for a class of equation systems that incorporate 
linear operators. A superior iteration method has been developed by deliberately shunning an attractive 
natural relationship that involves those linear operators in favor of concocting a clumsy artificial relationship 
that involves only function values. The point of proceeding in this way is that purely arithmetic operations 
with function values require no approximations^ whereas even quite elementary- looking operations involving 
operators may not be practically feasible without making use of potentially disastrous perturbation expan- 
sions. Indeed function manipulations can, on the contrary, be directed toward the goal of achieving iteration 
schemes that have continued fraction rather than perturbational character. 



Conclusion 

We have constructed a simple spherically-symmetric and energetically self-consistent Newtonian gravitostatic 
model which yields a simple and apparently very useful upper bound on the amount of self-gravitationally 
corrected energy which can be contained in a spherical region. That bound is just the "Planck force" (c^/G) 
times the radius r of the sphere, a relationship which, inter alia, implies that the Schwarzschild radius 
never lies in free space, making the Schwarzschild singularity physically unrealizable. This same model 
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can also be obtained by using the Birkhoff theorem and its corollary to extend the "Newtonian" form of 
the free-space Schwarzschild metric into the interior region of its self-gravitationally corrected spherically- 
symmetric static source. We have, moreover, been able to extend this Newtownian gravitostatic model to 
nonsymmetric energy distributions by means of an interesting nonperturbational continued-fraction iteration- 
solution method. That extension turns out to yield qualitatively the same upper bound (i.e., to within a factor 
of two) on the amount of self-gravitationally corrected static energy which can be contained in a spherical 
region. 

Such a "Planck force" times radius upper bound on a sphere's contained energy cuts off the mass deviation 
of an interacting quantum virtual particle at approximately the Planck mass, which in principle does away 
with the ultraviolet divergences that bedevil quantum particle scattering amplitude calculations. 

The "Planck force" times radius bound on a sphere's contained energy ought to be attained for contained 
quantum fields, which have infinite energy before self-gravitational correction, due to the combination of their 
unbounded frequency spectra and the quantum uncertainty principle. But only the universe itself, with its 
cosmological redshift, is actually capable of "containing" the arbitrarily high frequencies of a quantum field. 
Roughly estimating the radius of the universe as its age times the speed of light, and then dividing the 
"Planck force" times this radius by the corresponding spherical volume that has this radius, yields a rough 
averaged "dark energy" density estimate of about 1.7 joules per cubic kilometer, which is of the same order 
of magnitude as observational data. The same formula suggests that the early universe might have had an 
immensely greater "dark energy" density, perhaps as much as a Planck unit of energy density, which would 
be roughly 120 orders of magnitude times its present value. It is interesting that this seems to provide 
an automatic mechanism for the inflation of the early universe. The formula also suggests that the "dark 
energy" density will be decreasing toward zero as the universe expands, but that it won't decrease as rapidly 
as the density of ordinary matter will, which will increase the relative dominance of dark energy. 

Finally, there seems to be no reason why the presently evolved "dark energy" density should not share the 
small inhomogeneities which are so typical of the rest of the universe, such as the small peaks in the cosmic 
microwave background, and the galaxies, groups, filaments and voids in the distribution of luminous matter. 
If "dark energy" indeed has inhomogeneities, then might not those inhomogeneities themselves be the thing 
we call "dark matter" ? In spite of all the gravitational evidence for "dark matter" , there is apparently no 
observationally- known non- gravitational signal whatsoever for it. It would be a relief if something so elusive 
ultimately turned out to not have an independent existence. 
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